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ABSTRACT 

Globular clusters are an important test bed for Newtonian gravity in the weak- 
acceleration regime, which is vital to our understanding of the nature of the gravita- 
tional interaction. Recent claims have been made that the velocity dispersion profiles 
of globular clusters flatten out at large radii, despite an apparent paucity of dark mat- 
ter in such objects, indicating the need for a modification of gravitational theories. 
We continue our investigation of this claim, with the largest spectral samples ever 
obtained of 47 Tucanae and M55. Furthermore, this large sample allows for an accu- 
rate metallicity calibration based on the equivalent widths of the calcium triplet lines 
and K band magnitude of the Tip of the Red Giant Branch. Assuming an isother- 
mal distribution, the rotations of each cluster are also measured with both clusters 
exhibiting clear rotation signatures. The global velocity dispersions of NGC 121 and 
Kron 3, two globular clusters in the Small Magellanic Cloud, are also calculated. By 
applying a simple dynamical model to the velocity dispersion profiles of 47 Tuc and 
M55, we calculate their mass-to-light profiles, total masses and central velocity disper- 
sions. We find no statistically significant fiattening of the velocity dispersion at large 
radii for M55, and a marked increase in the profile of 47 Tuc for radii greater than 
approximately half the tidal radius. We interpret this increase as an evaporation sig- 
nature, indicating that 47 Tuc is undergoing, or has undergone, core-collapse, but find 
no requirement for dark matter or a modification of gravitational theories in cither 
cluster. 
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1 INTRODUCTION 

The nature of the gravitational interaction is one of the 
most important concepts in astrophysics, yet complete com- 
prehension of this interaction is still elusive. The so-called 
Pioneer and Fly-by anomalies, where spacecraft exhibit be- 
haviour that is unexpected from Newtonian and general rel- 
ativity gravitation theories, outline this lack of understand- 
ing (see l Anderson et al. Il2002l : lde DiegolbOOSl . and references 
therein), although these examples may have more mundane 
explanations. More importantly, it has been claimed that 
several globular clusters (GCs; u Centauri, M15, M30, M92 
and M107) exhibit a flattening of their velocity dispersion 
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profiles at radii R where r t is the tidal radius of the 

cluster (|Scarpa et aLl2003l . [2004al lbl). The authors claim that 
either dark matter (DM), or a modification of gravitational 
theory, is required to explain their results. 

Modified theories of gravity (MOG; see 
iDurrer fc Maartend 1200^ . for a review of modified gravity 
theories) and those of Newtonian dynamics (MOND; 
iMikroml Il983l l have been shown to solve some of the 
discrepancies. However, these are not universal theories 
and to-date have only been applied to specific instances 
(e.g. the Bullet C l uster and galaxy rotation cu rves; see 
lAngus et ai\ l2006l : ISanders fc NoordermeeJ |2007| . respec- 
tively). MOG theories diverge from Newtonian gravity in 
the high-acceleration regime and MOND diverges from 
Newton in the low-acceleration regime. Therefore, if either 
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of these theories were correct, the effect should be measur- 
able at the predicted accelerations. Independent of MOG 
or MOND theories, testing the gravitational interaction at 
low accelerations is essential to the overall understanding 
of gravity. 

Globular clusters are an ideal testing ground for weak- 
field gravitation because the accelerations experienced by 
stars at large radii are below the limit where DM, or a 
modified gravitation theory, is required to explain observa- 
tions in maiiy dyn amical systems (ao ~ 1.2 x 10~^'^ms~^; 
IScarpa et a/.ll2007l ') . Furthermore, they are thought to con- 
tain little, or no, dark matter - indicated by dynamica l 



models (|Phinnevlll993h . A^-body simulations (IMoor 



observations of GC tidal tails (jOdenkirchen etal. 



199e), 



2001) 



dyna mical and luminous masses of GCs ( Mandushev et alA 
Il99ll ) and th e lack of microlensing events from GC- mass 
dark haloes (iNavarro et al. I Il997l : llbata et al. 1 120021 '). al- 
though this is still under debate. GCs are also located at 
varying distances from the centre of the Galaxy, so that if 
all exhibit similar behaviour. Galactic infiuences cannot be 

the pri mary cause. 

In [l ane et al ] I2OO9I (hereafter IPaoer J ) we calculated 
the velocity dispersions and mass-to-light profiles of M22, 
M30, M53 and M68. Our conclusions were that there is no 
requirement for significant quantities of dark matter, or a 
modification of Newtonian gravity, to explain the kinemat- 
ics of any of these clusters. In the current paper we continue 
this investigation with the largest spectroscopic dataset to 
date of 47 Tucanae and M55. We begin by describing the 
data acquisition/reduction (Section [2]) and the membership 
selection for each cluster (Section I2.2|l . Our large samples 
allow us to determine a metallicity calibration based on the 
Tip of the Red Giant Branch (TRGB; Section [STT} , as weU 
as the rotations fSection l3.2|) . systemic velocities and veloc- 
ity dispersions ( Section 13. 3 |l . and mass-to-light profiles (Sec- 
tion 13. 4p - where we also place limits on the DM content 
of each cluster from their velocity dispersions and mass-to- 
light profiles. Finally, our concluding remarks are presented 
in Section |4l 



2 DATA ACQUISITION AND REDUCTION 

AAOmega, a double-beam, multi-object spectrograph on the 
3.9m Anglo- Australian Telescope (AAT) at Siding Spring 
Observatory in New South Wales, Australia, was employed 
to obtain the data for this survey. AAOmega is capable of 
obtaining spectra for 392 individual objects over a two de- 
gree field of view. We used the D1700 grating, which has 
been optimized for the Ca II infrared triplet region, centered 
on 8570A, with 30 sky fibres used for optimal sky subtrac- 
tion, and 5-8 fibres for guiding. The positional information 
for our targets was taken from the 2MASS Point Source 



Catalogue |skrutskie et a/.ll2006l ) which has an accuracy of 
^ 0.1", and we selected stars that matched the J — K colour 
and K magnitude range of the red giant branch (RGB) of 
each cluster. 

Our observations were performed over 7 nights on Au- 
gust 12-18 2006, and a further 8 nights on August 30 - 
September 6 2007. The mean seeing was ~ 1.5". Several 
fibre configurations were taken for each cluster with 3600- 
5400 second exposures giving a signal-to- noise of ~ 50 — 250. 



To minimize scattered light cross-talk between fibres, each 
field configuration was limited to stars in a 3 magnitude 
range. In total, 4670 and 7462 spectra were obtained in the 
47 Tuc and M55 regions, respectively. Flat field and arc lamp 
exposures were used to ensure accurate data reduction and 
wavelength calibration. The pointing accuracy of the AAT 
is ~ 0.3" and the fibres have a ~ 2" diameter. The offset 
due to the pointing uncertainty is azimuthally scrambled by 
the fibre, so has no effect on the zero point of the wave- 
length calibration. Data reduction was performed with the 
2df dr pipelin«[3, which is specifically designed for AAOmega 
data. The efficacy of the pipeline has been checked with a 
comparison of individual stellar spectra. 

Radial velocity and atmospheric parameters were ob- 
tained through an iterative process which tak es the best 
fits to synthetic spectra from the library by iMunari et al\ 
(2005) and cross-correlates this model with the observed 
spectra to calculate the radial velocity [a process very simi- 
lar to that used by the Radial Velocity Ex periment (RAVE; 
ISteinmetz et aZ.lbOOelJz witter et alll2008l 'l projec t]. We used 
the same spectral library as the RAVE studies; iKiss et "all 
(|2007i ) outline this process in detail. 

2.1 Radial Velocity and Uncertainty Estimates 

To be sure that we are not under/overestimating the un- 
certainties, and that our radial velocity measurements, and 
estimates of the random uncertainties, are reproducible, we 
observed a s ingle con figuration of M68 (as part of the data 
obtained for lPaperJ ). consisting of 317 spectra of the same 
stars on consecutive nights. 

Radial velocities from the data for a single night were 
estimated using two independent pipelines to test the effi- 
cacy of our own software. For this we compare d the outputs 
of our own pipeline (again see lKiss et al ]|2007l . for a detailed 
description) with that from the pipeline written specifically 
for the RAVE project. The average difference in radial ve- 
locities between the two pipelines is 0.3 ± 0.1 kms^^. 

To test that our pipeline reproduces reliable velocities, 
and associated uncertainties, these were extracted from the 
data from consecutive nights. Subtracting one from the other 
results in a distribution with a mean of — 0.33kms~^ and 
a « 0.97. Therefore 95% of the stars observed have veloci- 
ties within « 2kms~^ from one night to the next. This is 
within the systematics of the instrument (see Section I3.3|l 
and is comparable to the quoted uncertainties for individ- 
ual velocity estimates. Furthermore, because we used many 
fibre configurations for each cluster, it was necessary to test 
for systematic offsets between configurations. Therefore, we 
calculated the velocity dispersion of stars from four separate 
configurations within the same distance bin (13 stars were 
available from each configuration). The dispersions between 
configurations had a maximal difference of 0.3 kms~^, which 
is well within the uncertainties of the bin. 



2.2 Cluster Membership 

We selected cluster members using four parameters, namely 
the equivalent width of the calcium triplet lines, surface 
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gravity (log g), radial velocity and metallicity ([m/H]). Only 
stars that matched all criteria were judged to be members. 
A further cut of log g < 3.25 was applied to 47 Tuc to ensure 
as many Galactic stars were as possible were removed. This 
very probably removed some cluster members but was a nec- 
essary sacrifice to ensure our sample was as free of Galactic 
contaminants as possible. Figure [1] shows the selection crite- 
ria of 47 Tuc. The selections of M55 are not shown because 
the velocity of the cluster (~ 175kms~^) is far removed 
from the Galactic population (~ Okms^^) so the selections 
are very clean. 

It s hould be noted that for several clusters studied in 
IPaper j . a cutoff of TcS > 9000 K was necessary to remove 
hot horizontal branch (HB) stars with highly uncertain ra- 
dial velocities because the calcium triplet in very hot stars 
is replaced by hydrogen Paschen lines (and also have intrin- 
sic radial velocity variability, see Section |3.3|I . No cut was 
made on T^g for either of the current clusters because no 
stars with > 7050 K (for 47 Tuc) or Toff > 8500 K (for 
M55) remained after our selection process. Figure [2] shows 
the relative locations of the observed stars and highlights 
those found to be members. Several member stars in each 
cluster were found beyond the tidal radius; the implications 
of this are discussed in Section [3311] 

Figure |3] shows the colour-magnitude diagrams (CMDs) 
of the cluster members, with the extra-tidal stars as large 
points. These extra-tidal stars do not populate any specific 
region of the CMD, indicating that there is no systematic 
contributing to their selection as members. The HB, RGB 
bump and asymptotic giant branch (AGB) clump are all 
visible in 47 Tuc. A hint of the HB can be seen in M55 at 
~ 13 < J < 14, however, it is not well populated because, 
unlike 47 Tuc, M55 has a blue HB, and these stars are too 
hot to exhibit strong calcium triplet spectra. A total of 2241 
and 726 members were selected for 47 Tuc and M55, respec- 
tively. For 47 Tuc, 98.6% of our final sample used for analy- 
sis (2210 out of 2241 stars) fall within the 99.7% confidence 
level for cluster membership based on statistical analysis of 
each selection parameter. 



2.2.1 Sgr, NGC 121 and Kron 3 

M55 resides in front of the Southern tidal ta il of the Sagit- 
tarius dwarf galaxy (Sgr; llbata et al\ 1 19941 ') and 42 stars 
from that field were found to be part of Sgr (see Section 
13.11 for details). These were removed from the sample and 
analysed separately. Two GCs from the Small Magellanic 
Cloud (SMC) are present in the 47 Tuc field, namely NGC 
121 and Kron 3. When determining Kron 3 membership, 
we restricted our selection to a box centered on the clus- 
ter with a length of 204 " , the diameter of the cluster as 
quoted bv iBica fc Dutral (|2000l ) (Figure [J). For NGC 121 
we simply selected the clump of stars shown in Figure [1] 
because there were no stars outside ~ 150", despite the di- 
ameter of this cluster being about the same as that of Kron 3 
(IBica fc Dutrall2000D . We then overplotted the selections on 
the 2MASS region around 47 Tuc to ensure they fell in the 
same region of sky as the cluster and discarded those that 
did not (Figure [4]). A total of 10 and 11 stars were found to 
be members of NGC 121 and Kron 3, respectively. 
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Figure 4. Selected members of Kron 3 (large points below 47 Tuc: 
(RA,DEC)ai(6. 19, -72.79)) and NGC 121 (large points above 47 
Tuc: (RA,DEC)r;(6. 70,-71. 54)). 



3 RESULTS 
3.1 Metallicity 

Given the large sample sizes of both clusters we were able 
to perform an accurate metalli city calibration in a simila r 
way to that bv lCole et all (|2004l ') and lWarren fc Cold (|2009D . 
Our method uses the TRGB K m agnitude {Ktrgb ), instead 
of the HB used by ICole et al\ |2004i ') and 3Varren & Coli 
(|2009l '). Both methods are robust, but the TRGB is brighter 
than the HB, allowing our method to be used for more dis- 
tant clusters where the HB is not visible. In addition, only in 
the K band has a direct c alibration of the TRG been made 
with Hipparcos parallaxes (|Tabur et a/.i.2 009V The metallic- 
ity calibration was carried out in three steps. Firstly, Ktrgb 
was subtracted from each star and plotted against the equiv- 
alent width of the calcium triplet lines (see Figure [5]), giving 
a distance independent measure of the lum inosity. Ktrgb 
values were ta ken from IValeriti et al\ l|2004h (47 Tuc, M30 
M55 and M681 . iMarconi et al\ l| 19981 ) fSgr core) and 2MASS 
CMDs (Kron 3, NGC 121, M22 and M53; see below). Sev- 
eral stars from the M55 field overlapped the 47 Tuc region 
in this figure. Because the metallici ty of 47 Tuc is s imilar 
to that of the tidal tails of Sgr (e.g. IChou et aT]|2007f ). this 
overlap was expected; these stars belong to the Southern 
tidal tail of Sgr and are highlighted in Figure [5] All stars 
with Toff > 6000 K were removed from the M55 sample for 
the metallicity analysis because these hotter HB stars have 
the calcium triplet lines affected by hydrogen Paschen lines 
and, therefore, should not be used for metallicity analysis. 
This was not necessary for 47 Tuc since the selections were 
restricted to stars with logp < 3.25, because this cluster is 
closer to the Galaxy in logg, which removed all the hotter 
HB stars. 

The second step in the calibration process was to fit 
straight lines to the data (Figure (5)1 . Because our 47 Tuc 
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Figure 1. The selections made for 47 Tuc. Tlie boxes indicate the selections on each parameter. The lower left panel also shows the 
selections for Kron 3 (at ~ 2600" radius) and NGC 121 (at ~ 2100" radius). The selection box for Kron 3 is centered on the cluster and 
was restricted to the cluster diameter (204"). For NGC 121 there were no stars outside 150" despite the diameter being similar to that 
of Kron 3. 




AX (°) AX (°) 

Figure 2. The uncircled points are the stars which were observed and determined not to be cluster members. The members used for 
analysis, bas ed on the selection method described in the text, are circled points. The large circle is the tidal radius of the cluster from 
iHarrisI l|l996l V In each panel. North is up and East is to the left. 
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Figure 3. CMDs of selected members of 47 Tuc and M55 with extra-tidal stars shown as large points. Because these do not populate 
any particular part of the CMD, there is no systematic contributing to their selection as members. Note the AGB clump (A), horizontal 
branch (B) and RGB bump (C) in 47 Tuc. The HB of M55 at ~ 13 < J < 14 is sparsely populated because it has a blue HB whose stars 
are too hot for strong calcium triplet spectra. 



CO - 



CD 



C\2 - 





47 Tuc 

M55 ' 



4 



K 



K. 



RGB 



tip 



Figure 5. K — Ktrgb vs equivalent width of the calcium triplet 
lines. Crosses are M55 members, dots are 47 Tuc members and 
circled crosses are those stars determined to belong to the South- 
ern tidal tail of Sgr. These were analysed separately. The straight 
lines are the linear fits to the data once the Sgr members, and all 
stars having Tgff > 6000 K, were removed. 



sample contains many HB stars, and these were not ex- 
pected to exhibit the same relation between calcium triplet 
line widths and K — Ktrgb , these fits were also performed 
with the HB stars of 47 Tuc removed. No difference in the 
fits was found, so the HB stars are included in Figure (5] 



The slope of these lines are 0.47 and 0.42 for 47 Tuc and 
M 55, respectively ( c,f. 0.64 for V band analysis of the HB 
bvlCole et al\\2QQi a nd 0.47 for K band analysis of the HB 
IWarren fc Colell2009l .'l. 



Thirdly, by plotting [Fe/H] of the two clusters (from 
iHarrid Il996l ) vs Y,W - AX, where A is the gradient of 
the slope above and X K ~ Ktrgb, f or these two clus- 
ters we have a calibrator on [Fe/H] (see ICole et Zll2004l : 
IWarren fc Coli 2009, for a detailed discussion of this cali- 
bration methodology). TW — AX can then be calculated for 
any cluster and therefore [Fe/H]. Figure [6] shows [Fe/H] cal- 
culated by this method versus [Fe /H] from th e literature: for 



Sgr ( Chou ct al. 2 007,) Kron 3 (|Glatt et al . 2008b), NGC 
121 dGlatt et aLll2008al) as well as the fou r cluster s from 



IPaperJ, na mely M22 (iMonaco et all |2004|) M30 (|Harrisl 



1996), M53 (|Harrislll996h and M68 ((Harridil996l ) (TRGB 
values for these final four clusters were all measured from 
2MASS CMDs), showing this calibration for [Fe/H] to be 
robust. For clarity, the [Fe/H] values from this paper, and 
the literature, are also shown in Table [1] 

The large uncertainty for Sgr is due to the difficulty 
in determining an accurate measure of the TRGB for this 
object, because the Sgr tidal ta ils are close to the G alaxy 
in colour- magnitude space (e.g. iMarconi et fflZ.|[l998l V Fur- 
thermore, there is a metallicity gradient along the tails. To 
determine the K magnitude of the TRGB for the South- 
ern Sgr tail at the location covered in the current survey, 
we produced a CMD from our data (Figure (Tj and found 
Ktrgb ~ 11.15. The paucity of stars in this CMD intro- 
duces large uncertainties into the calculated [Fe/H] for this 
object (Figure [6]). 

Since no literature values of the Ktrgb are available 
for Kron 3, NGC 121, M22 or M53, we produced CMDs 
from 2MASS of stars within 2' of the cluster core (for M22 
and M53) and within 10" for Kron 3 and NGC 121. Several 
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[Fe/H] (this paper) 

Figure 6. Calculated [Fe/H] for M30 (square), M68 (star), M53 
(circle), M22 (triangle), NGC 121 (cross), Sgr (diamond), Kron 
3 (asterisk) versus literature values (see text). The large uncer- 
tainty for Sgr is due to the difficulty in calculating the TRGB for 
this object. No consensus on an uncertainty estimate has been 
reached, hence none is shown here from the literature. 



Table 1. [Fc/H] values from this paper and from the literature, 
in order of decreasing mctallic ity. The li terature values ar e taken 
from : Kron 3 (Glatt et al. 2od^}j_Ss£j[Cho^^O[iQQ3)jJ;JGC 
121 jGlatt et a l. 2008a), M22 (Monaco ct al."2004'), M53 l|Harrisl 
119961) . Meslllarris 1996) and M30 (Harris 1996). 



Cluster 


[Fc/H] (this paper) 


[Fe/H] (literature) 


Kron 3 


-1.05 ±0.1 


-1.08 ±0.12 


Sgr 


-1.4 ±0.5 


1.1 


NGC 121 


-1.5 ±0.1 


-1.46 ±0.1 


M22 


-1.78 ±0.15 


-1.68 ±0.15 


M53 


-1.99 ±0.1 


-1.99 ±0.08 


M68 


-2.06 ±0.15 


-2.06 ±0.11 


M30 


-2.16 ±0.15 


-2.13 ±0.13 



hundred stars in the CMDs for M22 and M53, and ~ 100 
for Kron 3 and NGC 121, meant that accurate values could 
be calculated. 



3.2 Rotation 

Assuming each cluster has an isothermal distribution, their 
rotations were measured by halving the cluster by position 
angle (PA) and calculating the mean radial velocity of each 
half. The two mean velocities were then subtracted. This was 
performed in steps of 10°, and the best-fitting sine function 
overplotted. The results are presented in Figure |51 

This method gives an amplitude that is twice the ac- 
tual measured rotation. Therefore, 47 Tuc exhibits rotation 
at 2.2±0.2 kms~^ with an approximate projected rotational 
axis along the line PA = 40° — 220°, and M55 shows rota- 
tion at a level of 0.25 ± 0.09 km s~^ and an approximate 
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Figure 7. CMD of Sgr stars extracted from the M55 field showing 
the TRGB a.t K ^ 11.15. 
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Figure 8. The rotation of each cluster calculated as the difference 
between the mean velocities on each side of the cluster along 
equal position angles, as described in the text. The best fit 
sine function is overplotted, and a typical error bar is represented 
in the lower left of each panel. 



Testing Newtonian Gravity with Globular Glusters 7 



47 Tuc 




6 

^ o 

^ CV! 



o 



10 



12 
J (mag) 



14 



16 



Figure 9. Velocity versus J magnitude of the members of 47 Tuc. 
The HB stars (B) appear to have the largest velocity dispersion 
of any stellar type. As per Figure |3] the AGB clump is labelled 
'A', the horizontal branch 'B' and the RGB bump 'C. 



axis of rotation along the line PA = 65° — 245°. Our rota- 
tio n measure for 47 Tuc cor responds to the value calculated 
by iM eylan fc Mayorl jldSdl ) at a radius of » 20 pc and by 
IStrug atskaya (lOsI) at « 36 pc. lSzekely et al\ ((2007, ) showed 
that M55 is rotating with a velocity of ~ 0. 5kms~^. This 
discre pancy can probably be attributed to ISzekely et al\ 
l|2007h having a sample size approximately half that of the 
current study. For both clusters, we corrected the individual 
stellar velocity data for the measured rotation before the 
velocity dispersions, and M/Ly profiles, were calculated. 




40 60 
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Figure 10. Velocities, and associated uncertainties, versus dis- 
tance of 47 Tuc. Complete sample (top panel) and HB stars only 
(bottom panel). Note that the HB stars are evenly distributed 
throughout the cluster. The horizontal line denotes the measured 
systemic velocity. 



3.3 Velocity Dispersions 

Figure [5] shows velocity versus J magnitude for the 47 Tuc 
data. It appears that the HB population (labelled 'B') has 
the greatest velocity dispersion of any stellar population in 
our sample. Stellar pulsations can alter the velocity disper- 
sion profile of a GO, if pulsating stars are present in suffi- 
cient numbers, so it is important to check for t his effect. 47 
Tuc only has one single known RR Lyrae star (|Bono et al\ 
I2OO8I ) so this will not affect the dispersion profile, however, 
it is natural to expect that the large number of HB stars in 
this sample to have an effect, since HB stars are known to 
pulsate when located close to the instability strip. To test 
this we calculated the velocity dispersion of the 47 Tuc sam- 
ple in J magnitude bins. The bin containing the HB stars 
(12.43 < J < 12.59) has a velocity dispersion of 7.5±0.3 and 
the overall dispersion is 7.5 ±0.6. Because these stars do not 
show an increase in dispersion compared with the complete 
sample, and, furthermore, are distributed evenly throughout 
the cluster (Figure [10)). we see no reason to exclude the HB 
stars from the velocity dispersion analysis. Our M55 sample 
contains very few HB stars (see Figure [Sjl so this effect is 
negligible in this cluster. 

We measured the systemic velocity of each cluster, us- 
ing a Markov Chain Monte C arlo (MCMC) maximum likeli- 
hood method jCregorvlbOOsI '). which takes into account the 



Table 2. Compari sons be t ween the systemic radial velocities of 
each cluster in the iHarrij lll996l ') catalogue and those from this 
survey. Velocities are in kms~^. 



Cluster 


Vr (Harris 1996) 


Vr (this paper) 


47 Tuc 


-18.7±0.2 


-16.85 ±0.16 


M55 


174.8±0.4 


177.37 ±0.13 



individual velocity uncertainties on the stars, providing the 
systemic velocity with associated uncertaintie s. A compari - 
son between our values and and those from the lHarri^ (|l996l l 
catalogue are shown in Tabled The velocities show system- 
atic differences of the order of 2 — 3 kms~^ between our mean 
values and t hose iuiHarr is (199§). These differences are sim- 
ilar to what lBalog et al\ (|2009l ) found for two open clusters, 
NGC 2451A and B, using the same instrument and analysis 
method and comparing data to velocities in the literature. 
Our interpretation is that there might be a systematic un- 
certainty in the zero-point of our velocity system. 

To determine the velocity dispersions of our samples, 
we binned the stars in annuli centered on the cluster cen- 
tre, ensuring approximately equal numbers of stars per bin 
(~ 50 for M55 and ~ 100 for 47 Tuc; see Table |3] for the 
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Table 3. Dimensions of the bins, and the number of stars in each, 
used in the velocity dispersion analysis. Only the inner bin edges 
are shown. Values given are parsecs from the centre of the cluster. 
The final bin extends to 80 pc for 47 Tuc and 35 pc for M55. 





47 Tuc 


M55 




Inner bin ed 


ge Stars per bin 


Inner bin edge 


Stars per bin 


0.000 


99 


0.00 


51 


2.131 


102 


1.41 


45 


3.048 


99 


1.98 


49 


4.031 


99 


2.62 


52 


4.834 


101 


3.14 


49 


5.698 


100 


3.76 


51 


6.416 


100 


4.29 


45 


7.210 


100 


4.71 
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bin dimensions). The MCMC method described above was 
then used to determine the dispersion in each bin, and the re- 
sulting velo city dispe r sion p rofiles were overplotted with the 
best-fitting |Plummeij (|l91lh model. The central velocity dis- 
persion and the scale radius (rs; containing half the cluster 
mass) were used for the fitting. The Plummer model allows 
for the calculation of the total cluster mass from the central 
velocity dispersion (cto) and via (sec Dcjonghc 1987 for a 
discussion of Plummer models and their application) : 



Mtot = 



SttG 



Twenty five extra -tidal stars w ere found to be members 
of 47 Tuc {ft ~ 56 pc: lHarrij|l996l ). The velocity dispersion 
of 47 Tuc shows a marked increase for i? > 28 pc (this is dis- 
cussed in Section I3.3.ip . Because of this increase in velocity 
dispersion in the outer regions, the outer two bins were ex- 
cluded from the Plummer model fitting to the dispersion 
profiles (Figure llip , since including them in the fit would 
have created an artificial increase in the fit over the entire 
cluster, and hence artificially altered the total mass, scale 
radius and central dispersion estimates. For M55, only 3 
extra-tidal stars were found to be members, and since there 
is no increeise in the dispersion profile, these are not affecting 
the profile in the outskirts of the cluster. The total masses, 
scale radii and central velocity dispersions are presented in 
Tabl e 131 Our mas s estimates agree well with other studies 
(e.g. lMevlanlll989l:[PiTor fc Mevlanlll993l : iMeziane fc ColinI 
ll996l : lKruiissen fc Mieskell2009l ), none of whom used Plum- 
mer models to calculate their estimates. 
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Figur e 11. Velocity d ispersion profiles of each cluster. The best 
fitting lPlummeJ ((m^) 

model is overplotted. Note that for 47 Tuc 
the outermost 2 bins are not included during the fitting process 
to ensure no artificial increase in the estimates of total mass, Ts 
or central dispersion. 



Table 4. Total masses, scale radii (r^) and central velocity dis- 
persions (o-q) for 47 Tuc and M55. 



Cluster 


Total mass (Mq) 


Ts (pc) 


o-Q (kms -*-) 


47 Tuc 


1.1 ±0.1 X lO'' 


7.8 ±0.9 


9.6 ±0.6 


M55 


1.4 ±0.5 X 10^ 


11.7±4.2 


2.7 ±0.5 



IScarpa et ai\ (|2003l . l2004al lbh showed an apparent flat- 
tening of the velocity dispersion profiles of five of the six 
GCs studied, indicating a significant DM component, or a 
modified theory of gravity, was required to explain their re- 
sults. MOG models generally differ from Newtonian gravity 
for large accelerations (e.g. in galaxy clusters or elliptical 
galaxies) but become Newtonian for intermediate accelera - 
tions [e.g. for solar system bodies; see lMoffat fc TothI (|2008l ) 
as an example]. MOND, however, becomes non-Newtonian 
for accelerations below about 1.2 x 10~^°ms~^ (|Milgroml 
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1 19831 ). approximately the regime where dark matter is in- 
voked to explain, for example, rotation curves of galaxies. 

Plummer models have the advantage of being monotoni- 
cally decreasing and, therefore, any flattening of the profiles 
would be discernible. Within the limits of the model, and 
the uncertainties in the data, neither velocity dispersion ex- 
hibits any flattening, although the proflle of 47 Tuc exhibits 
an increase in its dispersion. This could not be called "flat- 
tening" and a MOND/MOG model, or DM, is not required 
to explain this phenomenon, although a DM component may 
be one explanation (see Section [3.3.1[) . Therefore, we infer 
that neither dark matter, nor a modification to the current 
understanding of gravitation, are needed to explain the dy- 
n amics of e ither 47 Tuc or M55, corroborating earlier results 
in lPaper J for M22, M30, M53 and M 68, and similar conclu- 
sions drawn bv lSollima et all (|2009l ') for u) Centauri. 

Because we only sampled 10 and 11 stars from NGC 
121 and Kron 3 respectively, it was not possible to create a 
velocity dispersion profile for these objects. Instead, a single 
velocity dispersion value for the cluster was calculated, with 
NGC 121 having a velocity dispersion of 2.2 ± l.lkms"^ 
and Kron 3 of 1.8±0.9kms~\ 

3. 3. 1 Evaporation 

GCs are known to be tidally dest royed by their host galax y 
through tidal heating (e.g. Pal 5: lOdenkirchen et aZ.l20C)a ). 
The tidal stripping signature is evidenced by stars being 
stripped in two directions (the leading and trailing tidal 
streams) . Because the sampled field around 47 Tuc does not 
reach far outside the tidal radius (Figure [2}, it is not possi- 
ble to tell whether the extra-tidal stars exhibit a preferen- 
tial direction. Therefore, an inspection of a 20° x 20° region 
centered on 47 Tuc was performed, using 2MASS data se- 
lected on the RGB and HB colours/magnitudes of 47 Tuc. 
No evidence for extended tidal stru cture was observe d, in- 
dependently confirming the result bv lLeon et al\ l|200(t ) who 
found no convincing statistical evidence of tidal tails. 

Evaporation in GCs has been shown to occur due to 
internal two-body relaxa tion over the cluster lifetime (e.g. 
iMcLaughlin fc FallllioO^ ). particularly during core-collapse, 
or in post-core-collapsed clusters. Importantly, A'^-body sim- 
ulations have shown this evaporation exhibits a signature in 
the velocity dispersion profile, increasing the dispersion at 
rt/2, precisely the r egion where our 47 Tuc profile increases 
l|Drukier et aZ.ll2007l : for 47 Tuc, rt/2 ^ 28 pc). Furthermore, 
the extra-tidal velocity distribution is symmetric about the 
systemic velocity (Figure llOp implying that these stars are 
being accelerat ed in a symme t ric wa y into the outer regions 
of the cluster. iDrukier ef all (|2007l ) pointed out that this 
evaporation is exacerbated by the collapse of a GC core, 
leading to greater numbers of two-body interactions. This 
leads to greater numbers of stars accelerated to the outer 
regions of the cluster, and beyond, increasing evaporation. It 
is this evaporation scenario that we infer for 47 Tuc, from its 
velocity dispersion at large radii, adding to the growing body 
of evidence that 47 Tuc is in a core-collap s e, or post-core- 
collapse, phase (e.g. [ Cebhardt et al\ Il995l : iRobinson et al\ 
ll995l : lHoweU et aZ.lbOOOl ." and references therein). 

To confirm that our interpretation is valid, and that 
we are not simply describing a chance phenomenon, we re- 
binned the data to have 50 stars per bin. The velocity 



dispersion profile is unaltered by the different binning. Fur- 
thermore, we changed the bin boundaries to be sure that this 
increasing velocity dispersion is real, and, again, found no 
difference in the overall shape of the profile. Of course, alter- 
native explanations exist for the increase in dispersion. For 
example if GCs form in a similar fashion to Ultra Compact 
Dwarf galaxies, there may be a large quantity of DM in the 
outskirts of the cluster as discussed bv lBaumgardt fc Mieskd 
(2008). However, no evidence exists supporting GCs form- 
ing in this manner. A full MCMC analysis of the outermost, 
apparently evaporating, stars will be performed in a subse- 
quent paper. 



3.4 Mass-to-Light Profiles 

In dynamical systems such as elliptical and dwarf galax- 
ies, large quantities of dark matter are evidenced by high 
mass-to-light ratios (M/Ly ^1). DM causes higher stellar 
accelerations, leading to inherently higher maximal stellar 
velocities, and hence a larger velocity dispersion. Therefore, 
one method for determining whether a pressure-supported 
object like a GC is DM dominated is to measure its M/Ly 
from its velocity dispersion. For our M /Ly profiles we hav e 
used the surface brightness profiles bv lTrager et al\ l|l995l l. 
converted to solar luminosities per square parscc, and our 
density profiles were calculated using (|Deionghail987i '): 



p{r) = 



Mtot 



(r2 -hr2)2 



The M/Ly profiles, and associated M/Ly values, are shown 
in Figure 1121 Because of the uncertainty in core mass, the 
stated M/Ly values are the mean mass-to- light for R > rs- 
Our method for calculating M/Ly is preferable to the widely 
adopted method using the central mass and luminosity, due 
to the uncertainty in core mass. Despite its apparent supe- 
riority over other methods for measuring M/Ly, very few 
studies have adopte d it. T his technique has been used by 
I Cebhardt fc Fischer! (|l995l ) for 47 Tuc, but only for the in- 
ner 10' (< rs). Because the mass estimates are highly un- 
certain at small radii we do not claim to have any realistic 
data on the M/Ly for R < rs, so no comparison b e tween 
the current study and that by ICebhardt fc Fischer! (!l995! ) 
can be made. 

Neither cluster has M/Ly S> 1, therefore DM cannot 
dominate, although the larger M/Ly of 47 Tuc may indi- 
cate a small DM component. However, it has been shown 
that ultra-compact dwarf galaxies, which follow the same 
luminosity - velocity dispersion relati on as old GCs, show 
no evidence of DM for M/Ly < 5 (e.g. !Hasegan et aZ.![2005! : 
!Evstigneeva et ttZl!2007l ). Because of this, we see no require- 
ment for DM in either cluster. 



4 CONCLUSIONS 

Having the largest sample of spectra ever obtained for 47 
Tuc and M55, we were able to produce a very accurate cal- 
ibration of [Fe/H] based on the equivalent width of the cal- 
cium triplet lines and the K band m agnitude of the TR GB. 
This method is similar to that by !Cole et aZ] |20o3) and 
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Figure 12. Mass-to-light profiles of 47 Tuc and M55. The thick 
hnc is the calculated M/Ly and the thin lines are the uncertain- 
ties. The vertical line is r^, and the mean M/Ly is only calculated 
for R > Ts due to the uncertainty in core mass. Neither cluster 
has M/Ly 2> 1, indicating dark matter is not dominant. 



IWarren fc Coi3 (|2009D . except that we use the TRGB in- 
stead of the HB which means this method can be used for 
much more distant objects. 

We calculated the rotation of our clusters assuming 
them to be isothermal. The rotation of 47 Tuc is ~ 2.2 ± 0.2 
with an approximate projected rotational axis along the line 
PA = 40° - 220°, and M55 exhibits rotation at a level of 
0.25 ± 0.09 kms"^ and has an approximate axis of rotation 
along the hne PA = 65° - 245°. For 47 Tuc, the rotation 
amplitude is in good agreement with previous work (e.g. 
iMevlan fc Mavoij|l986t). The only previ ous study estimating 
the rotation of M55 (jSzekelv et al .|2007^ found a value about 
twice that of the current work. but lSzekelv et al\ (|2007l 'l had 
a sample size approximately half that of ours, which may ex- 
plain this discrepancy. 

Our calculated velocity dispersion profiles of 47 Tuc 
and M55 provide no evidence that either DM or a mod- 
ification of current gravitational theories are required to 
reconcil e their ki nematic properties, corroborating previous 
work in iPaper ll . T he dynam i cs of M55 are well described 
by a purely analvtic lPlummeil (jToiJ) model, which indicates 



that Newtonian gravity adequately describes its velocity dis- 
persions, and shows no breakdown of Newtonian gravity at 
ao ~ 1.2 X 10"^" ms~^, as has been claimed for other GCs. 
The internal dynamics of 47 Tuc (for R < rt /2) are also very 
well described by the Plummer model, however, the veloc- 
ity dispersion profile of 47 Tuc exhibits a large increase for 
R > rt/2, exactly the region where evaporation due to two- 
body interactions in the core should be observable, especially 
for GCs undergoing core-collapse or in a post-core-collapse 
state. We interpret this increase in velocity dispersion as 
evaporation, and hence that this cluster is either presently 
in a state of core-collapse, or is a post-core-collapse GC. 
This adds to the growing evidence that 47 Tuc is currently 



unde r going a dynamical ph a se change (e.g . Gebhardt et al\ 
1 19951 : Ikobinson eTai] Il995l : iHowell ei ad l2000l . and refer- 
ences therein). A full analysis of the outer regions of this 
apparently evaporating cluster will be performed in a sub- 
sequent paper. 

We used a Plummer model to determine the total mass, 
scale radius (r^), and the M/Ly profile for each cluster. We 
find that neither cluster has M/Ly S> 1, indicating that 
DM is not dominant. Within the uncertainties, our esti- 
mated cluster masses matc h those in the literature well, 
as do the M/Ly ratios ( e.g. lMevlan lll989l: iPrvor fc MevlanI 
1 19931 : iMeziane fc ColinI Il996l : iKruiisscn fc Mieskd |2009|). 
The mass-to- light profiles produced bv iGebhardt fc Fischerl 
(1995) cannot be compared to the current work because they 
sampled the inner 10' for which the mass is uncertain. Note 
that we consider using mass-to-light profiles is a more accu- 
rate method for calculating M/Ly than using the core mass 
and surface brightness, because of this uncertainty. 

While our results strongly indicate that the current un- 
derstanding of globular clusters being dark matter poor, and 
their dynamics explained by standard Newtonian gravity, 
more robust dynamical modelling is required for confirma- 
tion. 
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